Constraining the nuclear gluon distribution in eA processes at RHIC 
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A systematic determination of the gluon distribution is of fundamental interest in understanding 
the parton structure of nuclei and the QCD dynamics. Currently, the behavior of this distribution at 
QQ ' small X (high energy) is completely undefined. In this paper we analyze the possibility of constraining 

' the nuclear effects present in xg^ using the inclusive observables which would be measured in the 

, future electron-nucleus collider at RHIC. We demonstrate that the study of nuclear longitudinal 

^Nj ■ and charm structure functions allows to estimate the magnitude of shadowing and antishadowing 

; I ' effects in the nuclear gluon distribution. 

^ ■ PACS numbers: 12.38.-t, 24.85. +p, 25.30.-c 
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Since the early days of the parton model and of the first deep inelastic scattering (DIS) experiments, determining 
^ • the precise form of the gluon distribution of the nucleon has been a major goal of high energy hadron physics. Over the 
^ ' last 30 years enormous progress has been achieved. In particular, data from HERA allowed for a good determination 
, of the gluon density of the proton. A much harder task has been to determine the gluon distribution of nucleons 
bound in a nucleus, i.e., the nuclear gluon distribution {xg^{x,Q^)). In recent years several experiments have been 
', dedicated to high precision measurements of deep inelastic lepton scattering (DIS) off nuclei. Experiments at CERN 
J> ■ and Fermilab focus especially on the region of small values of the Bjorken variable x = Q'^ /2Mv, where = —<f 
' is the squared four- momentum transfer, v the energy transfer and M the nucleon mass. The data [l|, taken over a 
. wide kinematic range 10~^ < x < 0.1 and 0.05 Gel^^ < < 100 Gel^^, show a systematic reduction of the nuclear 
] structure function F^{x, Q'^)/A with respect to the free nucleon structure function F^{x,Q'^). This phenomenon is 
, known as nuclear shadowing effect and is associated to the modification of the target parton distributions so that 
' xq^{x,Q'^) < Axq^{x,Q^), as expected from a superposition of pp interactions (For a review see, e.g. 0,3). The 
modifications depend on the parton momentum fraction: for momentum fractions x < 0.1 (shadowing region) and 
0.3 < X < 0.7 (EMC region), a depletion is observed in the nuclear structure functions. These two regions are 
bridged by an enhancement known as antishadowing for 0.1 < a; < 0.3. The experimental data for the nuclear 
structure function determine the behavior of the nuclear quark distributions, while the behavior of the nuclear gluon 
distribution is indirectly determined using the momentum sum rule as a constraint and/or studying the logQ^ slope 
of the ratio F2"' / F2 Currently, the behavior oi xg^{x,Q'^) at small x (high energy) is completely uncertain as 
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d shown in Fig. [U where we present the ratio Rg = xg^ / {A.xg^), for A — 208, predicted by four different groups 
which realize a global analysis of the nuclear experimental data using the DGLAP evolution equations [5| in order to 
determine the parton densities in nuclei. In particular, the magnitude of shadowing and the presence or not of the 
antishadowing effect is completely undefined. 

In the last years the analysis of the nuclear effects in deep inelastic scattering (DIS) has been extensively discussed 
01 and motivated by the perspective that in a near future an experimental investigation of the nuclear shadowing 
at small x and » 1 GeV^ using eA scattering could be performed at Brookhaven National Laboratory (eRHIC). 
It is expected that measurements over the extended x and ranges, which would become possible at eRHIC, will 
give more information in order to discriminate between the distinct models of shadowing and the understanding of 
the QCD dynamics at small x. This collider is expected to have statistics high enough to allow for the determination 
of several inclusive and exclusive observables which are directly dependent on the behavior of the nuclear gluon 
distribution, as for example, the longitudinal and charm structure functions, the logarithmic slopes with respect to x 
and Q^, as well as the diffractive leptoproduction of vector mesons. In particular, the longitudinal structure function 
is expected to be measured for the first time in the kinematical regime of small x, since the electron - ion collider will 
be able to vary the energies of both the electron and ion beams. 

In this paper we study the behavior of the nuclear longitudinal structure function F^ and the charm structure func- 
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tion i^2 '^ and analyze the possibility to constrain the nuclear effects present in xg^ using these inclusive observables. 
We estimate the normalized ratios 

considering four distinct parameterizations for the nuclear gluon distributions and compare their behavior with those 
predicted for the ratio Rg = xg^ / A.xg^ . We analyze the similarity between these ratios and demonstrate that the 
experimental study of these observables allow to determine the magnitude of shadowing and antishadowing effects. 
We calculate these observables using the Altarelli-Martinelli equation [s^ and the boson-gluon fusion cross section 
0, respectively. In other words, we will restrict ourselves to the descriptions which use the DGLAP evolution 
equations Q to describe the behavior of the nuclear parton distributions and will assume the validity of the coUinear 
factorization. It is important to emphasize that the theoretical understandingof small- a; and large A regime of the 
QCD dynamics has progressed in recent years (For recent reviews see, e.g. [lOj), with the main prediction being a 
transition of the linear regime described by the DGLAP dynamics to a nonlinear regime where the physical process 
of parton recombination becomes important in the parton cascade and the evolution is given by a nonlinear evolution 
equation. One of the main motivations for the eRHIC experiment is the study of this new regime, denoted Color 
Glass Condensate (CGC) ^T^. As in Ref. 11 1 the inclusive observables at eRHIC were studied using a generalized 



saturation model, based on the CGC physics, the current study can be considered as complementary to that reference. 

Let us start presenting a brief review of the calculations of the longitudinal and charm structure functions. The 
longitudinal structure function in deep inelastic scattering is one of the observables from which the gluon distribution 
can be unfolded. Currently, there is a expectation for new experimental HERA data for taken with reduced 
proton energies, which will provide more direct access to the proton gluon distribution and shed light on the QCD 
dynamics at small- a; (See, e.g. Ref. [l^). Longitudinal photons have zero helicity and can exist only virtually. In 
the Quark-Parton Model (QPM), helicity conservation of the electromagnetic vertex yields the Callan-Gross relation, 
Fl = 0, for scattering on quarks with spin 1/2. This does not hold when the quarks acquire transverse momenta from 
QCD radiation. Instead, QCD yields the Altarelli-Martinelli equation|8j] 

Fl{x, Q') = x^ f %[\ F,iy, + 4 ^ e^(l - ^)yg{y, Q^)] , (2) 

-^T^ J X y "J q y 

expliciting the dependence of Fl on the strong coupling constant and the gluon density. At small x the second term 
with the gluon distribution is the dominant one. In Ref. ISj the authors have suggested that expression ([2]) can 
be reasonably approximated by Fl{x,Q^) « 0.3 ^^xg{2.5x,Q'^), which demonstrates the close relation between the 
longitudinal structure function and the gluon distribution. Therefore, we expect the longitudinal structure function 
to be sensitive to nuclear effects. In this paper we calculate F]^ using the Altarelli-Martinelli equation 

Let us now discuss charm production and its contribution to the structure function. In the last years, both the 
HI and ZEUS collaborations have measured the charm component F2 of the structure function at small x and have 
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found it to be a large (approximately 25%) fraction of the total [Tj|. This is in sharp contrast to what is found at 
large x, where typically F2/F2 ~ 0(10^^). This behavior is directly related to the growth of the gluon distribution at 
small-x. In order to estimate the charm contribution to the structure function we consider the formalism advocated 
in 15] where the charm quark is treated as a heavy quark and its contribution is given by fixed-order perturbation 
theory. This involves the computation of the boson-gluon fusion process. A cc pair can be created by boson-gluon 

fusion when the squared invariant mass of the hadronic final state is > 4m^. Since — ^IT^"* + Mf^, where 
Mjv is the nucleon mass, the charm production can occur well below the threshold, w 4to?, at small x. The 
charm contribution to the proton/nucleus structure function, in leading order (LO), is given by [9| 



1 ,:.c/ ^2 2n o 2"^(/^'^) dy .X 
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where a = 1 + and the factorization scale [i! is assumed ^'^ — 4m^. Ci 2 is the coefficient function given by 



+ ;3hl + 8z(l-z)-z(l-z)^]}, (4) 

where /3 = 1 — q^^^I^-^ is the velocity of one of the charm quarks in the boson-gluon center-of-mass frame. Therefore, 
in leading order, 0(0^), F2 is directly sensitive only to the gluon density via the well-known Bethe-Heitler process 
7*5 — > cc. The dominant uncertainty in the QCD calculations arises from the uncertainty in the charm quark mass. 
In this paper we assume = 1.5 GeV . 

Finally, let us briefly discuss the distinct parameterizations for the nuclear parton distributions (For details see the 
recent review Q). We will make use of the existing parameterizations of the nuclear parton distribution functions based 
on a global fit of the nuclear data using the DGLAP evolution equations. Currently there are four parameterizations, 
proposed by Eskola, Kolhinen and Salgado [Ts'l, by de Florian and Sassot ^19*1, by Hirai, S. Kumano and T. H. Nagai 
|2(j | and the very recent one proposed by K. J. Eskola, H. Paukkunen and C. A. Salgado ^21]. In what follows they 
will be called EKS, DS, HKN and EPS, respectively. The basic idea of these approaches is that the experimental 
results presenting nuclear shadowing effects can be described using the DGLAP evolution equations with adjusted 
initial parton distributions. Similarly to the global analyzes of parton distributions in the free proton, they determine 
the nuclear parton densities at a wide range of x and through their perturbative DGLAP evolution by using the 
available experimental data from I A DIS and -pA collisions as a constraint. As pointed out in Ref. different 
approaches differ in the form of the parameterizations at the initial scale, in the use of different sets of experimental 
data, in the order of the DGLAP evolution, in the different nucleon parton densities used in the analysis, in the 
treatment of isospin effects and in the use of sum rules as additional constraints for the evolution. For instance, the 
DS and HKN groups provide leading (LO) and next-to-leading order (NLO) parameterizations, while EKS and EPS 
perform only a LO QCD global analysis. There are noticeable differences between the HKN analysis results and the 
ones in Ref. (loj especially in the strange-quark and gluon modifications. These differences come from various sources. 
First, the analysed experimental data sets are slightly different. Second, the strange-quark distributions are created 
by the DGLAP evolution by assuming s(x) = at the initial scale, and the charm distributions are neglected in 
Ref. [1!||. These differences lead to discrepancies among the gluon modifications. In contrast to the EKS, DS and 
HKN parameterizations, the EPS one has included the RHIC data from [13] in the global fitting procedure. The main 
assumption is that these data can be understood with linear evolution. The inclusion of the high-pj^ hadron data from 
RHIC at forward rapidities provided important further constraints for the gluon shadowing region. By construction, 
these parameterizations describe the current experimental data. However, the resulting parton distribution sets are 
very distinct. In particular, the predictions of the different groups for Rg differ largely about the magnitude of the 
shadowing and the presence or not of the antishadowing. It is associated to the fact that the data included in the 
global analyses probe the quark distribution, while the gluon is constrained only by the evolution and the momentum 
sum rule. As shown in Fig. [TJ while the HKN and DS parameterizations predict a small value of shadowing, the EKS 
and EPS one predict a large amount, with the distinct predictions differing by a factor A at x — 10^^. Furthermore, 
while the DS parameterization does not predict antishadowing and EMC effects in the nuclear gluon distribution, 
these effects are present in the EKS and EPS parameterizations. In the particular case of the HKN parameterization, 
it predicts a steep growth of the ratio Rg in the region x > 10^^. It is important to emphasize that the magnitude of 
shadowing and antishadowing effects in the EKS and EPS parameterizations are directly related by the momentum 
sum rule. The large discrepancies between the predictions of the four parameterization for xg"^ in all kinematical x 
range imply a large uncertainty in the predictions for the observables which would be measured in pA/AA collisions 
at LHC, for instance. 
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FIG. 2: Ratios Rg, Rc and Rl for the four considered nuclear parameterizations and = 2.5 GeV^. 



As mentioned above it is well known that the inclusive observables Fl and are strongly dependent on the 
gluon distribution. Our goal is to quantify and determine the kinematical region where these observables directly 
determine the behavior of Rg. In order to obtain model independent conclusions we calculate Rl and Rc using the 
four parameterizations described above and compare with the corresponding predictions for Rg . As the small- a; region 
at eRHIC will be probed at small-Q^ we concentrate our analysis on two characteristic values of Q^: = 2.5 GeV^ 
and 10 Ge'V^. Moreover, we only consider A — 208, but similar conclusions are obtained for othes values of the atomic 
number. 

In Figs, m and [3] we present our results. Firstly, let us discuss the small-x region, x < 10~^, determined by 
shadowing effects. We observe that Rl practically coincides with Rg for all parameterizations and for the two values 
of considered. This suggests that shadowing effects can be easily constrained at eRHIC by measuring Fl. This 
conclusion is, to a good extent, model independent. On the other hand, the ratio Rc gives us an upper bound for 
the magnitude of the shadowing effects. For example, if it is found that Rc is equal to w 0.6 at a; = lO^"' and 
= 2.5 GeV^ the nuclear gluon distributions from DS and HKN parameterizations are very large and should be 
modified. At — 10 GeV^ the behavior of Rc is almost identical to Rg, which implies that by measuring F2 at this 
virtuality we can also constrain the shadowing effects. Considering now the kinematical range of a; > 10^^ we can 
analyze the correlation between the behavior of Rl and Rc and the antishadowing present or not in the nuclear gluon 
distribution. Similarly to observed at small values of x, the behavior of Rl is very close to the Rg one in the large-x 
range. In particular, the presence of antishadowing in xg^ directly implies an enhancement in F^. It is almost 10% 
smaller in magnitude that the enhancement predicted for xg^ by the EKS and EPS parameterizations. Inversely, if we 
assume the nonexistence of the antishadowing in the nuclear gluon distribution at a; < 10^^, as in the DS and HKN 
parameterizations, no enhancement will be present in F^ in this kinematical region. Therefore, it suggests that also 
the antishadowing effects can be easily constrained at eRHIC measuring Fl. On the other hand, in this kinematical 
range the behavior of Rc is distinct of Rg at a same x. However, we observe that the behavior of Rc at x — 10^^ is 
directly associated to Rg at x — 10"^. In other words, the antishadowing is shifted in Rc by approximately one order 
of magnitude in x. For example, the large growth of -Rg predicted by the HKN parameterization at a; > 10^^ shown in 
Fig. [T] implies the steep behavior of Rc at x > 10~^ observed in Fig. O A similar conclusion can be drawn from Fig. 
121 Consequently, by measuring F2 it is also possible to constrain the existence and magnitude of the antishadowing 
effects. 

Some comments are in order here. Firstly, it is important to emphasize that although we have calculated Fl and F2 
at leading order we expect that the behavior of the ratios Rl and Rc and, consequently, the main conclusions of this 
paper would not be strongly modified by the NLO corrections. Secondly, in our study we only have considered two 
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examples of inclusive observables which would be measured at eRHIC. As demonstrated in Ref. [l6| the study of the 
logarithmic slope of the nuclear structure function is another important quantity to probe the nuclear effects and the 
QCD dynamics at small-x. Furthermore, the exclusive production of vector mesons is an important complementary 
observable to determine the nuclear gluon distribution, since in this case the total cross section is proportional to 
the square of xg^ (See e.g. Ref. [13 )• Finally, we have disregarded the presence of non- linear effects in the QCD 
dynamics and used the current parameterizations based on the DGLAP dynamics, extrapolating them to lower values 
of X. Consequently, our results can be regarded as conservative and serve as a baseline. Deviations from this baseline 
may indicate the emergence of the saturation regime of QCD. 

Summarizing, our results indicate that the study of the inclusive observables Fl and F| in eA process at eRHIC is 
ideal to constrain the nuclear effects present in the nuclear gluon distribution, which, in turn, is a crucial ingredient to 
estimate the cross sections of the processes which will be studied in the future accelerators. Basically, we see that by 
measuring these observables we will have a direct access to the nuclear gluon distribution and allow to discriminate 
between the different parameterizations. We hope that this paper can motivate a more accurate determination of 
and F2 in the next years. 
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